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Abstract. In this paper the perturbation theory with the frequency of transition in atom as perturbation 
parameter is constructed. The estimation of the reminder term of series of this perturbation theory is 
given. With the help of this perturbation theory we have found an exact asymptotics of eigenvalues of 
complete hamiltonian in the limit of high quantum numbers. It is shown that the counter-rotating terms 
keep a leading term but absolutely change a second term of this asymptotic. 


1. Introduction. 


The Jaynes-Cummings model without rotating wave approximation (RWA) is the elementary model 
describing an interaction of atom with a field. But despite of this it can not be solved exactly. This 
model without the RWA was considered by different methods in works [1-7]. The hamiltonian of this 
model has the form 

H = Hq + S' V = ujq Cq uj ci g ci (a-|- a^ ), (1) 


where a and a'*' are the photon creation and annihilation operators, g is the coupling constant, a; and 
cvo are the frequencies of mode and atomic transition respectively, (Tq and cti are the 2x2 matrices 
of form 


^2 0^ 

, (^1 = 

1—1 

o 

1—1 

o 




It is well known that the RWA formulas for eigenvalues take into account only zero and hrst order of 
the perturbation theory on the coupling constant g. Therefore they are valid only at small relative 
coupling constant g/u and sufficiently small quantum index. More precisely, the validity of the RWA 
formulas for eigenvalues is defined by the condition g^/n/oj -C 1. In the case of exact resonance 
(a; = (Uo) the expression gy/n defines the splitting of eigenvalues. In optics g/uj <^1. Hence, unique 
opportunity to leave for limits of the RWA is the consideration of the highly exited states with 
suffisiently large quantum index n. That is the RWA loses force at sufficiently large energies of a 
field mode. How the eigenvalues of the hamiltonian (|lD and the splitting of them behave at the 
large quantum indexes? In the present paper we shall answer on this question by constructing the 
perturbation theory on the parameter ujq, which enters linearly in the hamiltonian (|^). We shall 
show that this perturbation theory well describes not only an eigenvalues at (Uo/o; -C 1 but also the 
highly laying eigenvalues at arbitrary g and ujq < kuj, where k = -\/3/(27r) ~ 0.23 . We also give an 
estimation of the reminder term of series and find two first terms of asymptotic of eigenvalues on 
quantum index. It is interesting that the second term of this asymptotic is qualitatively differed from 




the corresponding term in the RWA. This difference leads to the fact that the splitting of eigenvalnes 
vanishes at the large qnantnm nnmbers, unlike the RWA case, when the splitting infinitely increases. 

We remain open the question about the validity of the Jaynes-Cummings model itself at the limit of 
large average energy of a field mode, because then the manyphoton transitions between other levels 
become important. Nevertheless, we can hope, that the resonant levels give the basic contribution in 
atomic dynamics even for large average energy of a held mode. 


2. Perturbation theory on the parameter ujq. 


In work we have shown that the hamiltonian of model (|l|) can be represented in invariant subspaces 
by two Jacobi matrices of form 
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Let us present the operators Hi and H 2 as 
Hi = Aq + cuo Pi , H2 = Aq + cuo P2 , 


( 2 ) 


where Aq is the unbounded main operator without periodic modulation of main diagonal. Pi and 






P 2 are the diagonal projectors 
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(4) 


The operator Aq is the hamiltonian of shifted oscillator. It can be diagonalized with the help of 
Bogolubov’s transformation. Its eigenvalues and eigenvectors have the form 

OO 

-^0 l^m) l^m) ; ^0 TflOJ Q jU) , |^m) ^ ^ ^ |^n) 7 ^7n,n 7 (5) 

n=0 

where {|erj,)} is the basis of matrix representation (|^), are defined by Feynman-Schwinger’s 

formulas [H,0 


Pi”l = exp{-j72a,n J— (») LZ'-'igyj) 

V ml Vo;/ 

Here are generalized Chebyshev-Laguerre’s polynomials 

(n + s)! 
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It is easy to veryfy that the expression (j^) can be presented also in the form of contour integral 

/// fgY-'^ 1 / /I ^ 


= exp{-//2.7 Y Y ^ f - 1) exp ^ ^ ^ d. 

c 


o;^ X 
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where C is the circle of unit radius with the centre in the origin of coordinates of a complex plane 
X. This expression we will use further. 













From (Q) it follows that the operators Hi and H 2 depend linearly on ujq. If we know the solution of 
spectral problem for the operator Aq, we can build the perturbation theory on the parameter ujq. 
Let us hnd the matrix form of the operators Pi and P 2 in the basis of the operator Aq eigenvectors. 
Let U( 5 ') be the orthogonal transformation from \em) at \am) 

\am)=V{g)\em), {ek\\J{g)\e^) = Pt\ U^U = E (8) 


Using (@),(|) and (||), we have 

= (a^lPilaJ = felWP, U|0 = (ei|P«|0 = P« pW (9) 

n—odd 

PS = {at|P2|a™> = fe|WP2U|e„> = (ei|pW|e.„>= i: pW pW, (10) 

n—even 


where P*^^^ = Pi U and P^^^ = P 2 U are the transformed projectors. 

Let us consider for example the sum (P). Using the representation and summarizing on odd values 
n, we come to the formula 


Pk]l = exp{-^V^^} ^k\m\ 
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The contour integrals in this expression can be calculated consistently with the help of residues. As 
a result, we obtain the following expression for 
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Similarly, one can obtain and the expression for P^}^^ dehned by the sum (p!0D 
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Comparing ( 0 ) and (|I^) with (^, we have 
1 


p(i) _ X ( ^) pM('o„'i 

- 2 - 


p(2) _ r I I -*-1 p(™)^o„'\ 


1 


(-1)^ 


(13) 


2 


(14) 














or in the operator form 


P(L = U^(^) P, U((/) = ^ (E - B \]{2g)) 

(15) 

P(2) = U^(^) P 2 U(^) = ^ (E + B U(2^)) , 

(16) 


where B is the diagonal matrix with elements Bm,k = (~1)^ ^m,k ■ Let us note that the matrices B 
and U satisfy to the identity 

[B U]2 = E 

The formulaes ( or 0 .( 0 ) ) allow to write at once the approximated expressions for 

eigenvalues taking into account only zero and hrst orders of the perturbation theory on ujq. The 
hrst order correcton to an eigenvalues is dehned by diagonal elements of perturbation. Taking into 
account the formula for the eigenvalues of the operator Aq and the expressions (|T3|), (^), (j^) (at 
k = m), we obtain the following approximated formulaes for eigenvalues and X^ of the operators 
Hi and H 2 respectively 

~ 3ujo/2 + muj- g^/u) - ^ exp |- ^| 

~ 3u;o/2 + muj - g^/uj + ^ exp |- ^| Lm{4:g^/uj^) 

This formulaes was obtained in work with the help of a so-called ’’operator method”. 

Let us consider now the constructed perturbation theory series in more detail. We shall show that 
the two hrst term of this series give an exact asymptotic of an eigenvalues X^ and X^^ at large 
quantum index m. 


3. Asymptotic of eigenvalues. 


Let us consider, for example, an eigenvalues A(^^ of the operator H 2 . The proof of the appropriate 
formulas for A(,(^ is completely similarly. In what follows for brevity we shall omit the top index (2) 
at eigenvalues and write Am instead of A(^\ The perturbation theory series for exact eigenvalue Am 
has the form 

00 

= ( 17 ) 

k=0 


General expression for A(^\ in case when an operator depends linearly on the perturbation parametr 
and an eigenvalues are not degenerate (here, due to the simplicity of Jacobi matrix spectrum), has 
the form []^ 
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Here, we have omitted as well as above the top index (2) at the perturbation operator and used 
the formula (|) for the unperturbed eigenvalues A^A 


We have found already that 

AJ,”> = mw + h)o - g^/oj 

^ ^ Wi)/2 + , m ^ CO 


( 20 ) 
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Here, we have used the asymptotic of Chebyshev-Laguerre’s polynomials (see, for example, H)- 
Let us consider the second order correction A(^^ which is dehned by 


a£> = 
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k^m 


(0) HO) 


A};i' - A 


According to ([I^) and (^0|), this expression can be presented in the form 


,(2) _ d V lPt\‘2g)? 

^ 4n) Ip m-k 

k^m 

The behaviour of this expression as m —> cx) is dehned by the behaviour of sum 
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k^m 


m — k 


( 22 ) 


Let us show that ^ 0 as m —> cx). For this purpose, let us transform (|22|) to the form 

00 ^ 

, V 

. 


n=l 
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where the transformation matrix Cm,n is dehned as follows 


r* — 


- [^i+m(2fi')]^ , n<m 
-[^i+m(2^)]^ n>m 


(23) 


The condition ^ 0 follows from ^ 0, if and only if the transformation Cm,n satishes to the 












following conditions ( jl^, Theorem 4) 


1. ^ \Cm,n\ < H, where H does not depend of m 

n 

2. lim Cjnn = 0 , for arbitrary n 

m^oo ’ 


( 24 ) 


Then is the regular transformation. Let us prove the hrst condition. Taking into account (p^), 
we have 


Elc'™,nl = E[[pi"l(2^)]'-[pi+U2(7)]2]+ E 
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.n=0 
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But since the values P^\2g) are the matrix elements of the orthogonal transformation \J{2g), the 
sum in square brackets is equel to unit identically. The diagonal matrix element P^\2g), due to 
(H), equals to 

Pln\‘^9) = 1 “ , 

and due to the asymptotic of Chebyshev-Laguerre’s polynomials 
Therefore, we have 

E I I < 1 ~ 0 as m ^ oo 

n 

And hence, the condition 1 in (p4D is fulhlled. 

Let us check now the validity of the second condition in p^). For this purpose, due to p^), it is nec¬ 
essary to consider the diagonal asymptotic of the non-diagonal matrix elements of the transformation 
\J{2g). From (|]), we have 
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Using the asymptotic of generalized Chebyshev-Laguerre’s polynomials JH 
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we obtain 
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Due to the symmetry of perturbation matrix Pk,m, we obtain at once the same asymptotic and for 

PSni^g) 


PSn{‘^g) - TJ7 , OO, for arbitrary n 


According to (p3D, it follows that the condition 2 is also fulhlled. Therefore due to the above theorem, 
tm and hence tend to zero as m —>• cxo : 

















Let us consider the third order correction to the eigenvalue 
that A^) is dehned by expression 


From 


and ([T^) it follows 
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Using (W), we have 
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Let us consider the hrst term in square brackets. Let us aplay Cauchy’s unequality to the sum on j 
in this composed 
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Due to the orthogonality of the transformation \J{2g), as well as above, we have 

'ji^rn = \/l — [Pj;^\2g)Y < 1, for arbitrary i and m , 

and therefore 
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Using the theorem on regular transformation, just as it was made for one can show that 
as m —>■ oo. 


Thus for the hrst term in we have the unequality 
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Let us apply once again Cauchy’s unequality to the sum on i in the right side of this unequality 
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It follows that 
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and the unequality (P6D takes the form 
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Since 'n j ~ 1.81 > 1 and P^\2g) 0 as m —>■ cx) that there exists such mg that 

1+ (-!)"■ Pir’(2s) <;r/V3, m>m„ 
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Despite of that \P^\2g)\ < 1 for arbitrary m, we could not prove that the unequality (1^) is valid 
for arbitrary m. Taking into account that a^n < 1 and (^), we obtain from (PS]) 
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Here, number 6 is the number of components in the sum (^) for k = 3. 


Let us note that since am —>■ 0 as m —>• cx), we could write instead of (|30| ) more strong unequality, 
following from (pSl) . But we shall write just the unequality (|30|), following from (^), because just 
this way can be used and for higher orders of a perturbation theory. 


Let us consider the k-th order correction Using Cauchy’s unequality and condition (^) one can 
show as well as above that the absolute value of each term in the sum (|^ is bounded by expression 
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where is the number of terms in the sum (^), i.e. the number of solutions of the equation 

Ui + ... + Ufc = fc — 1, Uj > 0 
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(33) 


Let US show the estimation (pT]), for example, on the typical term, entering in (^) at fc = 4 
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Using ([I^) , we have 
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At last, using (|^) and already obtained estimation (^), we obtain 
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that is the estimation (pT]) for k = A. 

With the help of the unequality (|3^) we can estimate the reminder term of the series (pAl) 
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From (p^ it follows that asymptoticaly 
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It follows that the series in right part of (^) converges at ujq < a;-\/3/(27r). We can estimate it as 
follows. Since W < we have 
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Taking into account that cXm —>■ 0 as m —a cx) and using (1^) ! (^) ) P^ ) and (E3), we obtain the 
following asymptotic of eigenvalues 


Am = mo; + 3ujq/2 — /oj + o(l), m —>■ oo 


(36) 


Since the formulas (P), ®) differ only by sign, it easy to see that the same asymptotic takes place 
and for eigenvalues A^^ of the operator Hi. Thus, we have proved the following result 


Theorem 


If can < a;\/3/(27r), then the eigenvalues A^^ and A^^ of the operators Hi and H 2 have the asymp¬ 
totic and the reminder term of perturbation theory series have the estimation (P). 
























4. Conclusion. 


In physical applications the main role plays not eigenvalues itself but a difference of neighbouring 
eigenvalues, determining in the resonant case Uq = u the splitting of originally degenerate levels 

~ ^2m+2 ~ ^2m.+l ) m = 0, 1, 2, . . . 

From (|36[) it follows directly that 

—>• a;, m —> cx) 


It is in the sharp contradiction with the RWA. In the RWA an eigenvalues appropriate, for example, 
to A(^) (in resonant case a; = a;o), are dehned by the expression 

Xjn = 00 (2m + 2) ± V2m + 1, m = 0,1, 2,... 


Therefore in the RWA the splitting grows as \/2m. 

This change of splitting undoubtedly should change the time dynamics of quantum amplitudes, 
especially, when the average energy of a held mode is sufficiently large. 

We have proved the asymptotic formula (^61) only at the condition coq < oo^j (27r). But the numerical 
calculations shows that it is valid and for cuq > a;\/3/(27r). 
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